Higgs Self-Coupling as a Probe of 
Electroweak Phase Transition 



Andrew Noble and Mcixim Perelstein 



Institute for High Energy Phenomenology 
QQ Newman Laboratory of Elementary Particle Physics, 

O Cornell University, Ithaca, NY 14853, USA 

^ anobleOphysics . Cornell . edu.maximOlepp . Cornell . edu 

1^ Abstract 

^ We argue that, within a broad class of extensions of the Standard Model, there is 

a tight corellation between the dynamics of the electroweak phase transition and the 
cubic self-coupling of the Higgs boson: Models which exhibit a strong first-order EWPT 
predict a large deviation of the Higgs self-coupling from the Standard Model prediction, 
as long as no accidental cancellations occur. Order-one deviations are typical. This 
shift would be observable at the Large Hadron Collider if the proposed luminosity 
or energy upgrades are realized, as well as at a future electron-positron collider such 
as the proposed International Linear Collider. These measurements would provide a 
laboratory test of the dynamics of the electroweak phase transition. 
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1 Introduction 



About 10^^" sec after the Big Bang, the universe underwent a phase transition in which the 
electroweak gauge symmetry was reduced from SU{2)l x U{1)y to t/(l)em- In the Standard 
Model (SM) and many of its extensions, the electroweak phase transition (EWPT) is under- 
stood in terms of the dynamics of a fundamental scalar field, the Higgs h: the vacuum h = 
preferred at high temperatures becomes unstable at a certain "critical temperature" and 
the Higgs develops a vacuum expectation value (vev) pQ. The nature of this transition is 
a subject of considerable interest. For example, is the EWPT a violent, out-of-equilibrium 
event accompanied by massive entropy production (first-order transition) or a gentle, quasi- 
adiabatic process (second-order transition)? Answering this question would have far-ranging 
implications: for example, electroweak baryogenesis, which theoretically is arguably the most 
attractive of the mechanisms proposed to explain the observed matter-antimatter asymme- 
try, can only occur if the EWPT is strongly first-order [2]. It is very difficult to obtain direct 
information on the EWPT from astrophysical observations, due to the long and complex his- 
tory of the Universe since the transition. (Note, however, that gravitational waves produced 
in a first-order transition might be observable at LISA [3j.) Luckily, the dynamics of the 
phase transition can be deduced theoretically, once the fundamental Lagrangian governing 
the physics at the electroweak scale is known. For example, it is well known that if there 
is no physics beyond the SM relevant at that scale, the transition is never first-order for 
the Higgs boson mass obeying the current experimental lower bound, 114.4 GeV. Of course, 
the gauge hierarchy problem and other arguments strongly suggest that physics beyond the 
Standard Model (BSM) should appear at the electroweak scale, and the nature of the EWPT 
depends on the (as yet unknown) details of that physics. Observing all of the weak-scale 
BSM particles and measuring their masses and interactions (in particular with the Higgs) 
would provide sufficient information to determine the order of the transition. In practice, 
however, this will be a very difficult task for the upcoming Large Hadron Collider (LHC), 
and even for a more distant next-generation e~^e~ collider such as the proposed Interna- 
tional Linear Collider (ILC), especially in models where some of the new states are gauge 
singlets. It would therefore be useful to identify simpler observables that are correlated with 
the order of the transition. In this paper, we propose that the value of the Higgs boson 
cubic self-coupling is one such observable. Specifically, we argue that models leading to a 
first-order phase transition predict a large (typically order-one) deviation of the Higgs cubic 
coupling from its SM value. We show that this is the case in three toy models, which capture 
three basic mechanisms for obtaining a first-order EWPT: loop-induced corrections to the 
Higgs effective potential, nonrenormalizable operators in the potential, and mixing between 
the weak-doublet Higgs and singlet scalars at the weak scale. The only exceptions to this 
statement that were observed in our analysis are due to accidental cancellations between 
large corrections to the Higgs cubic of different origins. Thus, while non-observation of a 
large deviation of the Higgs cubic from the SM prediction would not completely rule out a 
first-order EWPT, it would make it very unlikely. 

The correlation between a strong first-order EWPT and a large deviation of the Higgs self- 
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coupling from the SM value was already noted, in the specific contexts of low-cutoff models [1] 
and the two-Higgs doublet model [5]. Here, we will argue that this corellation is in fact much 
more general, and applies (with a small caveat requiring no accidental cancellations) in a 
very broad class of models. 



2 Renormalizable Models with Unmixed Higgs 

In this paper, we will not restrict ourselves to a specific extension of the SM. We will instead 
study a series of toy models, which capture the important features of realistic BSM theories 
exhibiting strong first-order EWPT. The toy models analyzed in this section satisfy two 
basic conditions. First, the only scalar field which changes its value during the EWPT is 
the neutral component h of the SM Higgs doublet^ H = {h + iG^)/y/2). With this 
assumption, the finite-temperature effective potential Ve^ih^T) completely specifies the dy- 
namics of the EWPT. Second, at tree level, the Higgs sector is weakly coupled and can be 
adequately described by a renormalizable Lagrangian, i.e. no irrelevant terms are important. 
Many theoretically motivated models satisfy these conditions, including the Minimal Super- 
symmetric Standard Model (MSSM) ^ for a broad range of parameters corresponding to 
the so-called Higgs decoupling region. Little Higgs models [7j such as the Littlest Higgs with 
T-parity, and others. Interesting models which violate either one of these conditions and 
exhibit first-order EWPT also exist; we will study such examples in sections [3] and |4} 
With the above conditions, the tree-level Higgs potential has the form 

V\h) = -^h^ + \h\ (1) 

while the masses of the SM and BSM particles are given by 

Ml{h)=Ml + a,h\ (2) 

At the one-loop level, the Higgs effective potential, including finite-temperature corrections, 
can be determined completely once the values of Mq^ and for every particle in the theory 
are known. In practice, particles with <C 1 (weak coupling to the Higgs) can be neglected. 
In the SM, these parameters are given by 

i = {t,W,Z,h,G), 
Ml = (0,0,0,-/i2,-/i'), 

- - ' ^,3A,A). (3) 
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^The "SM Higgs doublet" is defined as a field whose tree-level couplings to all SM states are given by 
the SM values, with no large corrections from new physics. For example, no significant tree-level mixing of 
h with gauge-singlet scalars is allowed. 
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The zero-temperature (Coleman- Weinberg) correction to the potential ([ij), in the MS regu- 
larization scheme, has the form 

i 

where A is the renormalization scale and Qi and Fi are the multiplicity and fermion number of 
the field i. In the SM, in the basis of Eq. Q , ^ = (12, 6, 3, 1, 3) and F = (-1, 1, 1, 1, 1). The 
constants Ci are regularization scheme-dependent and are physically irrelevant, since their 
only effect is an additive renormalization of the parameters /i^ and A and the cosmological 
constant. Rescaling A has the same effect; we fix A = 174 GeV and interpret /i^ and A as 
the values of the couplings at that scale. The parameters are subject to the constraint 

^ {V\h = v) + V^{h = v)) =0 (5) 
where v = 246 GeV. In addition, the Higgs mass, 

ml=^A^'(h = v) + V^{h = v)) , (6) 

must satisfy the experimental lower bound, 

nih > 114.4 GeV. (7) 

The potential defined in this way suffers from an infrared singularity that appears at the 
minimum of the tree-level potential, v^j.^^ = /i^/A, and is due to the massless (Goldstone) 
bosons G at that point [8]. While the potential itself remains finite {\imx-*o loga;^ = 0), the 



cubic coupling defined in (16) diverges at that point. Resummation of the infrared logarithms 
is required to correctly evaluate A3 in the neighborhood of that point. This resummation 
is conveniently achieved by switching from the zero-momentum renormalization scheme to 
the on-shell scheme [8]. Technically, this is equivalent to adding an extra term to Vq, and 
interpreting the second derivative at the minimum of the corrected potential as the physical 
(on-shell) Higgs mass. The additional term eliminates the infrared divergence. We will 
employ the on-shell scheme throughout this paper. 

The one-loop finite-temperature potential has the form [9] 

VriKT) = E is / - exp(-,8\/F + Mf{h))] 



-E 





9iT 
27r2 



dkk' log[l + exp(-^ Jfc2 + M^{h))] 



where f3 = 1/T. If the model contains bosons with masses small compared to T, multi-loop 
infrared-divergent contributions need to be included to correctly analyze the dynamics of h. 
These contributions can be resummed and written in the form of the "ring potential" [TO] : 



Vr{h, T) = 5^ ^ Tr [M^ih) - (Mfih) + 0,(0))^/^] , (9) 
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where the sum runs over all (light) bosonic degrees of freedom, and nj(0) is the zero- 
momentum polarization tensor. In the SM, 



n.(o) = uaio) = (A/ + 1,- + Ia? + Ia) , 

nGB(O) = ^TMiag(/,/,/,^'^) , 



(10) 



and 
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These expressions are valid at T ^ Mw', for smaller T, the ring terms are unimportant and 
can be dropped. The full potentiaj^is Ves{h; T) = + Vq + Vr + K- A useful "decoupling" 
feature of this potential is that at any T, the contributions from particles with Mi{h) ^ T 
are exponentially suppressed. Thus, in BSM theories, only weak-scale states need to be 
included to analyze the EWPT. 

At high temperature, T 3> Moi for all states, the potential has the simple form 



If the quantity in brackets is positive, (h) = is a stable minimum, and the electroweak 
symmetry is restored. This is the case, for example, in the SM and the MSSM, where > 
for all states. We will restrict our analysis to such models, since otherwise the EWPT does 
not occur. A first-order transition occurs if Ves{h; T) develops a second, symmetry-breaking 
local minimum, (h) = vt{T) ^ 0, while the high-temperature minimum (/;,) = is still clas- 
sically stable. Then, once the symmetry-breaking vacuum becomes energetically preferred, 
bubbles of the symmetry-breaking vacuum will start rapidly nucleating, growing, and coa- 
lescing, eventually encompassing the entire Universe and completing the phase transition. 
The critical temperature at which the phase transition occurs can be found from the 



■^For some values of h, the potential V(h,T^ develops an imaginary part, due to the quantum instability 
of a classical state with spatially uniform value of h [TT\ This effect does not affect the dynamics of a 
first-order transition, as long as the rate of this instability is small compared to the electroweak scale. This 
condition is always satisfied in our analysis. 

•^Note that the actual transition temperature may be lower if the bubble nucleation rate F is strongly 
suppressed compared to its natural value, F ~ -y. Phase transition does not occur until Y ^ R ^ /M-pi. 
Given the large hierarchy v <^ -Mpi, we will not include this possibility in our analysis. 




(12) 



conditionE] 



Kff((/i) =0;T,) = V,^{vt{T,)-T,). 



(13) 
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The "strength" of the first-order transition (more precisely, the degree of deviation from 
quasiadiabatic evolution) can be characterized by the dimensionless parameter 

£ = ^. (14) 

For example, it is well known that successful electroweak baryogenesis requires ^ ^ 1. 

Due to the complicated form of the potential, the problem of finding Tc and vt{Tc) in 
a given model is in general not tractable analytically. We have developed a numerical code 
which is used to evaluate these quantities in all examples studied in this paper. 

In the SM, the EWPT is always second-order once the bound ([T]) is satisfied. Below, we 
analyze the situation in a few simple extensions of the SM. 

2.1 Single BSM Scalar 

Consider a toy model in which a single real scalar state is added to the SM, with the potential 

V =^M^S^ + C\H\^S^. (15) 

This is the most general renormalizable potential which respects a symmetry S —S. 
(Without this symmetry, tree-level mixing between h and 5* would generally occur. Such 
models will be discussed below.) Note that while we use the SUSY-inspired notation ("^ for 



the coefficient in Eq. (15), this coefficient can in general be negative, as long as S does not 
become tachyonic until h > v and the lifetime of the vacuum at {h = v, S = 0) is larger than 
the present age of the universe. (The runaway direction can be stabilized by the addition of 
a small S"^ term, which does not affect our analysis.) 

The scalar contribution to the one-loop effective potential is given by Eqs. (|8|, ([o]), 
where M{S) = Mq + ('^h'^, lls{0) = |T^C^ and gs = 1- In addition, the ring terms due to 
the h and G fields have to be modified by adding the scalar contribution to the polarization 
tensors, Allfi{0) = AIIg{0) = j^T^(^. (We neglect additional contributions to polarization 
tensors due to S self-coupling or coupling to other BSM singlets.) In order to identify 
the points with a first-order EWPT, we scanned the parameters of the modelj^ Mq = 
. . . 800 GeV, = — 4 . . . 8. We also scanned the parameters /i and A, making sure that the 
constraints ^ and ([T]) are satisfied; the scanned region corresponds to rrih = 114: A . . . 244.4 
GeV. (We did not find any points with higher rrih and a strong first-order EWPT.) For every 
point exhibiting a first-order EWPT, we compute the one-loop value of the physical Higgs 
boson cubic self-coupling, 

_ 1 £iV\h) + V^{h)) 
A3=Q U=.- (16) 



■*We include the models with negative values of in the scan in order to demonstrate that they exhibit 
the same correlation between the Higgs cubic coupling and the ^ parameter as the positive-C^ models. We 
imposed the constraint that S should not be tachyonic in the range h = . . . 300 GeV, but did not analyze 
the lifetime of the conventional EWSB vacuum for these models, so some of the negative-^^ points may be 
phenomenogically unacceptable. A detailed analysis of this issue is beyond the scope of this paper, and its 
outcome would not affect the main point of our argument in any way. 
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The results are presented in Fig. [TJ In the left panel, we plot the parameters ^ and A3 
(normalized to the SM one-loop value) for the points in our scan; in the right panel, we plot 
ruh and A3. It is clear from the plots that the Higgs self-coupling is significantly larger than 
the SM value for all points with a first-order EWPT in this model. If m/j is known from 
experiment (as it surely will be, by the time the experiments become sensitive to the Higgs 
self-coupling), the correlation between A3 and ^ becomes tighter, as illustrated in Fig. [2] 

Experimental measurement of A3 requires collecting a significant sample of events with 
two Higgs bosons in the final state. At the LHC, this coupling can be determined with 
20-30 % precision for the Higgs mass between 160 and 180 GeV, for integrated luminosity of 
3000 fb^^, which can be achieved if the luminosity upgrade to 10^^ cm~^s~^ is realized [12] . 
If an energy upgrade to a 200 TeV VHLC is realized, the coupling can be measured with 
8 — 25% precision for 150 GeV < nih < 200 GeV [13]. For lower Higgs masses, the VLHC 
could make a similarly precise measurement of A3 using rare decays [13] . Also, the 500-GeV 
ILC is expected to measure A3 with 20% accuracy for nih < 140 GeV and a luminosity of 
1 ab~^ [I5]. Our analysis shows that if a strong first-order EWPT occurs due to an extra 
scalar contribution at one loop, and the Higgs is sufficiently light, these experiments will 
discover a deviation of A3 from the SM prediction. Higher center-of-mass energy would be 
required if rrih is larger; note, however, that strong first-order EWPT in this setup requires 
rrih ^ 250 GeV, so a modest increase in ^/s would be sufficient to cover all of the interesting 
parameter space. 

In order to better understand this result, consider our model with /x = Mq = 0. In this 
case, the one-loop Higgs potential at zero temperature has the form 

where A can be expressed in terms of the coupling constants of the model, and the con- 
tributions from Higgs loops have been neglected. If A > (which occurs for sufficiently 
large ("^ in our model), this potential has a stable symmetry-breaking minimum at h = 
A exp(— 87r^A/y4 — 1/4), and a shallow local maximum at h = 0. This is of course just a 
realization of the purely radiative EWSB via dimensional transmutation (DT) a la Coleman 
and Weinberg [16]. If now a small positive Higgs mass term /i^ > is added, the potential 
around </> = is dominated by this mass term, and the maximum at the origin is turned into 
a minimum. At the same time, if /i^ is sufficiently small, the symmetry-breaking minimum 
is preserved, resulting in a potential with a "bump" between the two minima. A bumpy 
zero-temperature potential essentially guarantees that the EWPT will be of the first order, 
and most of the points with strong first-order EWPT in our scan (the blue points in Fig. Q 



have this feature. The Higgs self-coupling for the pure DT potential, Eq. (17), is 66% larger 
than the (tree-level) SM valu^for the same v and rrih: 

DT^5m| ^ 5 sM,tree .^g. 

^ Qv 3 ^ ^ ^ 



^The same relation holds in any model with approximate conformal symmetry in the Higgs sector broken 
by nearly marginal operators |17j . 
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Figure 1: SM with a single extra scalar. Models with a "bumpy" zero-temperature Higgs 
potential are shown in blue, and those without the bump in red. (Left panel) The strength 



of the first-order EWPT ^, defined in Eq. (14), vs. Higgs cubic self-coupling. (Right panel) 



Higgs cubic self-coupling vs. Higgs mass for points exhibiting a strong first-order EWPT, 
^ > 1. In both plots, the Higgs self-coupling is normalized to the one- loop SM expectation 
for the same rrih. 



which is roughly at the center of the range for the models with strongly first-order EWPT in 
our scan. Even without the bump (red points in Fig. [T]), the models with ^ > 1 have zero-T 
potentials whose characteristic feature is a flatter maximum aX h = compared to the SM. 
This shape is close to the DT potential, again resulting in large deviations of A3 from the 
SM prediction. 

The connection with the pure DT case points to a potential concern about our model: 
In the DT case, the physical Higgs mass is given by rrih = A^/^v/(47r), so that large values 
of A are required to satisfy the experimental lower bound on rrih- Correspondingly, in our 
scans, fairly large values of ("^ (between 3 and 10) are required to obtain phenomenologically 
consistent points with a strong first-order EWPT. For ("^ ~ 10, one might worry about 
the validity of perturbation theory that was used in our analysis. Even for more moderate 

= 3 ... 5, two-loop corrections might be numerically significant. This is the case, for 
example, in the MSSM, where the stops effectively contribute with = ^/l2\f ^ 3.5 (see 



section 2.4). On the other hand, the main outcome of our analysis is not the precise values 



of A3 for each point in the scan, but the fact that the deviation from the SM in this quantity 
is large for models with strong first-order EWPT. It seems very unlikely that the higher- loop 
corrections would conspire to cancel the large one-loop correction to A3 and bring it back 
precisely to the SM value, so our main qualitative conclusion should continue to hold even 
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Figure 2: Same as left panel of Fig. [T} but with the Higgs mass fixed at 160 GeV. 



at strong coupling. 

It should be noted that it is possible to have a large deviation of A3 from the SM prediction 
without a strongly first-order transition. We found such points in our scan. These points 
have values of mo above, but not too far from, the Higgs vev v, and large values of a. 
The zero-temperature effective potential receives a correction of order a{v/mo)^, giving a 
substantial correction to A3; the correction to the finite-temperature potential decouples 
exponentially, as exp(— mo/T), and is not sufficient to yield a first-order transition. Thus, 
while an observation of a large correction to A3 would keep the door open to the possibility 
of a strongly first-order phase transition scenario, it would not prove it. 



Consider the SM with a single additional Weyl fermion x, with a Majorana mass term 



where dots denote higher-order nonrenormalizable terms. The mass has the form ([2]) with 
a = Mo/A, and the formalism presented above applies, with = 2. Adding a single 
fermion also modifies the ring contributions: An;,(0) = AnG(O) = ^T^a, where a is the 
same coefficient appearing in Eq. ([2]), and g = 2(A) for a Weyl (Dirac) fermion. We neglect 
additional ring contributions that would arise if x coupled to other BSM states. In this 
model, we performed the scan with mh = 114.4 ... 414.4 GeV, Mq = . . . 800 GeV, and 
a = —4 ... 8. For a > 0, A3 is always suppressed with respect to the SM value, while for 



2.2 Single BSM Fermion 




(19) 
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^ A3 A3 vs rrih for ^>1 
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Figure 3: SM with a single extra scalar-fermion pair. Models with a "bumpy" zero- 
temperature Higgs potential are shown in blue, and those without the bump in red. (Left 
panel) The strength of the first-order EWPT ^, defined in Eq. (14), vs. Higgs cubic self- 
coupling. (Right panel) Higgs cubic self-coupling vs. Higgs mass for points exhibiting a 
strong first-order EWPT, > 1. In both plots, the Higgs self-coupling is normalized to the 
one- loop SM expectation for the same rrih- 



a < 0, A3 is enhanced. (In the T = effective potential, an a < fermion behaves like an 
a > boson, and vice versa.) But in neither case do we find a strong > 1) first-order 
EWPT. 

Since the only inputs from the model into the EWPT analysis are the coefficients in the 
mass formula ([2]), the same analysis applies in a model with a single Dirac BSM fermion 
(with Mq = if the fermion is chiral, and possibly Mq 7^ if it is vector-like). The only 
change is the new multiplicity factor, = 4. This does not affect the conclusions. 

2.3 A Single BSM Scalar-Fermion Pair 

Any realistic extension of the SM has multiple bosonic and fermionic states, whose contribu- 
tions to the Higgs effective potential must be added. In this subsection, we will explore how 
the interplay between bosonic and fermionic contributions affects the connection between 
the Higgs cubic coupling and the strength of a first-order EWPT. To address this, consider 
a toy model with a single fermion (for concreteness, we assume it to be a vector-like Dirac 
fermion), and four identical scalars. The number of scalars is chosen so that the number of 
bosonic and fermionic degrees of freedom is the same, as would be the case in supersymmet- 
ric models. We searched for points with strong first-order EWPT in the following parameter 
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space: a, = = -1 . . . 4, M(f = . . . 600 GeV, = . . . 600 GeV, ruh = 114.4 . . . 294.4 
GeV. (We did not find any points with a strong first-order EWPT for higher m/j.) The 
results are shown in Fig. |3} In most of the parameter space, the correlation between ^ and 
A3, observed above in the purely scalar model, persists, and models with strong first-order 
EWPT predict observable enhancements of A3 with respect to the SM value. There is also 
a smaller region where an observable suppression of A3 is predicted; this is due to the dom- 
inance of the fermion loops, which contribute to A3 with a negative sign, for < Mq . 
There is also, however, an exceptional parameter region, in which A3 is close to the SM 
value. This region has the following origin. In the one-loop zero-temperature potential, 
fermion and scalar contributions cancel if Mq = M* and as = a^, since S and \E' form an 
exact supermultiplet, and A3 has its SM value. The scalar and fermion contributions to the 
finite-temperature part of the potential, however, do not cancel in this limit. In particular, 
if Qg is large, a large cubic term can be generated at high temperature, which can lead to 
a first-order transition. Fermion loops do not contribute to this term. This explanation 
highlights the exceptional nature of the region with SM-like A3: the relations = 
and as = a<i- cannot be enforced by a symmetry other than SUSY. Given that the S and 
fields have significant couplings to the Higgs sector, and that SUSY must be broken in that 
sector, these relations are not radiatively stable. For example, renormalization group flow 
will naturally induce a mass splitting between S and \E' of order 

=\M\S)-M\n -Tt-,Sml log (20) 

167r^ 

where Smh is the Higgs-higgsino mass splitting (typically of order 100 GeV or above) and 
A > 1 TeV is the messenger scale where logarithmic divergences are cut off. Given the large 
value of in the region where the first-order transition is possible, we expect SM to be 
naturally large, of order a few hundred GeV. Imposing that the mass splitting be at least 
100 GeV eliminates the region with first-order EWPT and no deviation in A3. We find that 
the minimal deviation in the cases 6M = (100, 200, 300, 500) GeV, is 7%, 10%, 17% and 
20%, respectively. We conclude that in the absence of accidental degeneracies in the mass 
spectrum, the A3 — ^ corellation works well in this model. 

2.4 Multiple BSM States 

So far, we have observed the correlation between A3 and ^ in simple toy models. Will it 
persist in realistic BSM theories with more complicated spectra and interactions? While 
it is clearly impossible to analyze all possible cases, it seems very likely that the answer is 
positive. Several arguments point in this direction. 

First, note that in a model with real (or A^/2 complex) scalar fields with identical 
masses and couplings, the zero-temperture potential (and therefore A3) is the same as that 
for an "effective" single real scalar with 

Mis = VnmI = v^C, 9es = 1. (21) 
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The same scaling applies for the scalar contribution to the finite-temperature potential (in- 
cluding rings), except one has to also rescale T N^^^T. If the scalar contribution dom- 
inates, this implies that the critical temperature in the A^-scalar model is reduced by N^^^ 
compared to the effective single-scalar model, and 0c is unchanged. We checked numeri- 
cally that this scaling is accurate to within 20% even after the SM contributions are added. 
Thanks to this scaling, the approximate correlation between ^ and A3 is preserved for > 1, 
and the models with > 1 and strong first-order EWPT predict larger minimal enhance- 
ment in A3 compared to the = 1 case. The same approximate scaling laws apply to models 
with N identical species of fermions as well. Since the gauge charges of the added states 
have no effect on the Higgs potential in the one-loop approximation, we conclude that the 
results of our toy-model analyses in fact apply to any model where a single multiplet (or a 
single supermultiplet) dominates the BSM contribution to the Higgs potential. The MSSM 
with light stops, and models with multiple gauge-singlet scalars as in Ref. [IH], are in this 
category. 

Even if multiple BSM states with different masses and couplings contribute to the Higgs 
potential, the quahtative picture should remain the same. In fact, we checked that a toy 
model with two scalars, whose masses and Higgs couplings are varied independently, produces 
results very similar to the single-scalar case. This is reasonable, since the similarity between 
the properties of the models with strong first-order EWPT and the DT model, noticed in 
section 2.1[ should persist independently of the precise origin of the large radative corrections 
to the Higgs potential. 

Based on this evidence, we believe that significant deviation of the Higgs cubic self- 
coupling from the SM value is a generic prediction of weakly-coupled models of EWSB with 
strong first-order EWPT, provided that no accidental cancellations occur. A measurement 
of this coupling can serve as a diagnostic tool to glean information about the nature of 
the EWPT. This would be particularly valuable in models with extra gauge singlets, whose 
masses and couplings to the Higgs are difficult to measure directly, especially if the decay 
h SS is kinematically forbidden. 



3 A Model with Nonrenormalizable Higgs Interactions 

The effect of new physics at the scale A ^ f on the Higgs potential can be parametrized by 
adding nonrenormahzable operators, suppressed by powers of A, to the potential. Grojean, 
Servant and Wells (GSW) pointed out that in the presence of such terms, strong first-order 
EWPT may occur [U |19l EO]. Following GSW, consider a model with a tree-level Higgs 
potential of the form 

V\H) = fi^H^H + X{H^H)^ + ^{H^Hf. (22) 

It is reasonable to use this effective theory to describe the EWPT if the electroweak symmetry- 
breaking vev is small compared to A. The strong first-order EWPT occurs for positive //^ 
and negative A. With these choices, it can be shown that the theory is applicable only if 
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Figure 4: SM with an operator. Models with a "bumpy" zero-temperature Higgs potential 
are shown in blue, and those without the bump in red. (Left panel) The strength of the 
first-order EWPT ^, defined in Eq. (14), vs. Higgs cubic self-coupling. (Right panel) Higgs 



cubic self-coupling vs. Higgs mass for points exhibiting a strong first-order EWPT, ^ > 1. 
In both plots, the Higgs self-coupling is normalized to the one-loop SM expectation for the 
same fUh- 



|A| ^ 1. In general, this requires fine-tuning: for example, if the dimension-6 operator arises 
due to exchanges of a heavy weak-singlet scalar, one can show that |A| ~ 1 will also be 
generated when the scalar is integrated out. This should be cancelled by a bare term in 
the lagrangian. However, given the relatively low values of A interesting for our analysis 
(500 GeV-1 TeV range), the required fine-tuning is numerically quite modest 

We analyzed the dynamics of the EWPT in the GSW model, following the general frame- 
work discussed in Section |2} We performed a scan of the following parameters: rrih = 
114.4 . . . 214.4 GeV, A = 100 . . . 1000 GeV. (We did not find any points with strong first- 
order EWPT for higher A and rrih-) Our results for the order of the transition are roughly 
consistentj^ with Refs. [H [20] . We also computed the Higgs cubic self-coupling (at the one- 
loop order) for each point in the scan. The results are shown in Fig. |4j In the region where 
there is a strong first-order EWPT, this model also predicts a large, order-one enhancement 

^Thc appropriate suppression of the quartic coupling can occur naturally in composite Higgs models [5T]. 

^Quantitative differences between our analysis and those of Refs. [H [50] arise due to different approxi- 
mations made in each analysis. Our analysis includes the full SM contributions to Vq, Vt and Vr, listed 
in Section [2] which were not included in Ref. [4 . On the other hand, we do not include the corrections 
to the one-loop potential and ring terms induced by the dimension-6 operator, and do not account for the 
possibility of delayed nucleation. Both effects were included in Ref. f20|. These differences do not affect the 
argument of this paper. 
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of A3 from its SM value, which will be easily observable at the ILC. This is easy to understand 
analytically: at tree level, we have 



^GSW 



1 + 



2v' 



(23) 



Since strong first-order EWPT can only occur for A in the 500 GeV-1 TeV range, the 



correction term is large for these points. (The one-loop corrections to Eq. (23) are at most 
about 10% in our scan.) In addition, as Fig. |4] indicates, there is a clear and strong correlation 
between this coupling and the strength of the EWPT. 



4 A Model with Tree-Level Higgs-Singlet Mixing 

In all examples studied above, the SM-like Higgs field h was the only field that acquired 
an expectation value during the EWPT. If other scalar fields are present, they may change 
their value in the same transition, so that the order parameter for the EWPT is effectively 
a linear combination of h and other fields. Of particular interest are models with extra 
gauge-singlet scalars, such as the NMSSM, certain Little Higgs models, and others. At low 
temperatures, both singlet and Higgs vevs are generically non-zero. At high temperatures, 
the Higgs vev is driven to zero, and although the singlet vev typically remains non-zero, 
the EWPT involves changes in both vevs. Since the zero-temperature singlet potential is 
not restricted to even-degree polynomials by gauge invariance, it may contain cubic terms, 
which can naturally produce large "bumps" and lead to strong first-order EWPT. We will 
consider a simple model of this kind in this section. Our main goal here is to illustrate that 
the connection between the strong first-order EWPT and observable deviations of the Higgs 
self-coupling from the SM value persists in this model. 

Many studies of models of this type have appeared in the literature, mostly in the context 
of the NMSSM and the nMSSM [22]. We will follow the conventions of a recent analysis by 
Profumo, Ramsey-Musolf, and Shaughnessy (PRS) [23], which applies more generally. (See 
also Ref. [25 •) Consider a model with a single real scalar S added to the SM. The most 
general renormalizable scalar potential has the form 

V\H, S) = -fi^H^H + A {H^H)' + '^H^HS + ^H^HS^ + ^S^ + ^ |^4^ ^24) 

where H = {H^, H^). We assume that the model has a stable vacuum at (H^) = v/\/2, {S) = 
s. (The conditions for this are listed in Ref. [23].) The physical Higgs fields are obtained by 
expanding around this minimum, H'^ = (h + v)/\/2 and S = s + x, and diagonalizing the 
mass matrix by a rotation 

hi = sin 9 s + cos 9 h , 

/12 = cos9 s ~ sin9 h . (25) 
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^ vs A3 for 




Figure 5: SM with a single extra scalar that mixes with the Higgs. Models with a "bumpy" 
zero-temperature Higgs potential are shown in blue, and those without the bump in red. 
(Left panel) The strength of the first-order EWPT ^, defined in Eq. (14), vs. Higgs cubic 
self-coupling. (Right panel) Higgs cubic self-coupling vs. Higgs mass for points exhibiting a 
strong first-order EWPT, ^ > 1. In both plots, the Higgs self-coupling is normalized to the 
tree- level SM expectation for the same rrih. 



The mixing angle is defined such that |cos6'| > With this convention, hi is more 

doublet-like and h2 is more singlet-like. For positive (negative) 6, the eigenmass m2 is 
greater than (less than) mi. 

We focus on the state hi, which for simplicity we will refer to simply as "the Higgs", and 
study the deviation of its cubic self-coupling from the SM prediction for a Higgs of the same 
mass. At tree level, 



AT! 

463 sin^ e + I264X sin^ 9 ] . (26) 



— ^ ( 12f A cos'^ 9 + 3ai cos^ 9 sin 9 + 6a2(x cos 9 + v sin 9) cos 9 sin 9 + 
6mf V 



Unlike the models studied in section [2| no large couplings are required to achieve first-order 
EWPT in this setup. We restrict our analysis to the points where all couplings are weak. 



and ignore the loop corrections to the relation (26). 

It is straighforward to compute the finite-temperature one-loop potential in this model, 
and analyze the dynamics of the EWPT numerically. We used the high-temperature approx- 
imation for y*, expanding it up to linear order in T. This approximation greatly speeds up 
the numerical analysis. We performed a partial scan of the parameter space of the model, set- 
ting a2 = and scanning oi = -150 ... - 50 GeV, 63 = -150 ... - 50 GeV, A = 0.15 . . . 0.25, 
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f VS A3 



Figure 6: The strength of the first-order EWPT C, vs. Higgs cubic self-couphng, in the model 



of Eq. (24). Only the mixing parameter, ai, is varied. The Higgs self-coupling is normalized 



to the one-loop SM expectation for the same irih- 



X = 150 . . . 250 GeV, and 64 = 0.5. This parameter region includes a large set of points 
with strong first-order EWPT, and lies roughly within the range consistent with precision 
electroweak constraint^ determined in [23j. The results of this scan are illustrated in Fig. [s] 
The models with strong first-order EWPT predict a suppression of the Higgs cubic self- 
coupling relative to the SM value, by a factor of 2 or more. The sign of the effect, however, 
is not unique in this model: for example, we found points with 02 7^ which predict a strong 
first-order EWPT and an enhanced Higgs self-coupling. The important point is that large, 
order-one deviations in A3 from the SM are typical. This is not surprising, since both these 
deviations and the possibility of strong first-order EWPT are due to same new tree-level 
interactions. 

The corellation between the size of this effect and the strength of the EWPT is less clear 
than in the examples of sections [2] and [3] This is simply due to a large number of parameters 
that affect A3, but have only a marginal effect on ^. If those parameters are fixed, and, for 
example, only the Higgs-singlet mixing parameter ai is varied, a clear corellation between A3 
and ^ emerges, as seen in Fig. |6} (In this figure we fixed 63 = —80 GeV, A = 0.20, x = 150 
GeV, 64 = 0.5, a2 = 0, and scanned Oi = —130 ... — 105 GeV.) 

Unlike the models of sections |2] and [3| where the non-standard Higgs cubic coupling 

^Note, however, that other new weak-scale states can contribute to precision electroweak observables and 
modify the fits, without changing the dynamics of the EWPT. A SM-like Higgs in the 300 — 800 GeV mass 
range is allowed if new physics produces a positive contribution to the T parameter. This is the case, for 
example, in the LHT model [25] , 
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might provide the first experimental hint of the existence of extra singlet scalars at the weak 
scale, the model studied here will provide a variety of deviations in the Higgs production 
cross section and branching ratios from the SM due to the singlet admixture ^3] • Many such 
deviations will probably be observed before the Higgs self-coupling measurement is available, 
especially if A3 is suppressed. Nevertheless, it is interesting that the general pattern of large 
deviations in A3 in all models exhibiting a strong first-order EWPT persists in this class of 
models. 



5 Conclusions 

In this paper, we investigated the link between the dynamics of the EWPT and the deviation 
of the cubic Higgs boson self-coupling in extensions of the SM. We studied this link in a series 
of toy models. Our toy models implement three mechanisms by which the weak-scale new 
physics could turn the second-order adiabatic transition predicted in the pure SM into a first- 
order transition. In the first set of models, the dynamics is determined by large quantum 
corrections to the Higgs potential, both zero- and finite-temperature, from loops of new 
particles (possibly, but not necessarily, gauge singlets) with strong coupling to the Higgs. 
In the second model, the first-order transition is due to the presence of a nonrenormalizable 
operator in the Higgs potential. In the third model, the crucial new feature is the tree- 
level mixing between the Higgs and an extra scalar singlet. In all cases, we found that 
models that exhibit a strong first-order EWPT predict large deviations of the Higgs cubic 
coupling from the Standard Model value. Order-one deviations are typical. The models of 
the first and second classes typically predict an enhancement in the Higgs cubic coupling, 
while in the models of the third class both signs are equally likely. The only cases where the 
cubic coupling does not deviate significantly from the SM prediction are due to accidental 
cancellations among unrelated contributions to this coupling. Thus, we conclude that a 
measurement of the Higgs cubic coupling consistent with the SM prediction would strongly 
disfavor (although not completely rule out) a first-order electroweak phase transition. A 
discovery of a substantial deviation of this coupling from the SM would keep the possibility 
of a first-order transition open, but would not conclusively prove it, as discussed at the end 
of section [2m 

While a measurement of the Higgs cubic coupling is very difficult at the LHC, a 20%-level 
measurement appears to be within the reach of the proposed luminosity or energy upgrades 
of the LHC, as well as the ILC. Based on the analyis, we conclude that such a measurement 
could serve as a sensitive probe of the dynamics of the electroweak phase tranition. If no 
deviation from the SM is observed, it would be very unlikely that the transition is first- 
order. Note that at least some of the sources of the first-order transition are very difficult 
to exclude by any other means: for example, a gauge-singlet scalar that does not mix with 



the Higgs at tree level but has a strong quadratic coupling, as in the model of section |2.1 
is essentially impossible to detect directly as long as the decay h SS is kinematically 
forbidden. Therefore, even with a large body of knowledge about the BSM particles and 
their couplings from the LHC and ILC experiments, the question of the order of the phase 
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transition would be very difficult to settle definitively. A measurement of the Higgs cubic 
coupling would be crucial in addressing this issue. 

It is exciting that the next generation of collider experiments will likely contribute to our 
understanding of cosmology. One example is the possibility that the dark matter particle will 
be discovered, and its propserties measured, at the LHC and the ILC. This connection has 
attracted much attention in the literature, and has been analyzed both in specific models of 
new physics [26] , and in a model-independent framework pT] . Understanding the dynamics 
of the electroweak phase transition is another open question in the early universe cosmology 
where the LHC and ILC experiments will have an important role to play. In this paper, 
we studied one simple experimentally observable quantity that is tightly correlated with 
the order of the transition. It is also possible that a combination of mass and coupling 
measurements could be used to completely reconstruct the Higgs potential, providing a 
much more detailed picture of the transition. It would be interesting to understand what 
quantities would need to be measured, and how precise the measurements would need to be. 
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